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EIGENVALUES FOR MAXWELL’S EQUATIONS WITH 
DISSIPATIVE BOUNDARY CONDITIONS 

FERRUCCIO COLOMBINI, VESSELIN PETKOV t AND JEFFREY RAUCH 

Abstract. Let V{t) = t > 0, be the semigroup generated by Maxwell’s 

equations in an exterior domain C with dissipative boundary condition 
Etan — 7 ( 3 ;)(i' A Btan) = 0 , 7 ( 0 :) > 0 , S F = We prove that if 7 ( 0 ;) is 
nowhere equal to 1, then for every 0 < e <C 1 and every £ N the eigenvalues 
of Gh lie in the region As U7?.]v, where As = {2 G C : | Re z\ < Cs(| Im 2 | 2 + 

1), Re 2 < 0}, TZ]^ = {2 G C : I Im 2 | < C]v(| Re 2 | + 1)“^, Re 2 < 0}. 


1. Introduction 


Suppose that A C {a; G : |a;| < a} is an open connected domain and 

il := \ K is an open connected domain with C°° smooth boundary F. Consider 

the boundary problem 

dtE = cm\B, (9(5 = —curl A in x fl, 

Etan - Btan) = 0 On K+ X F, (1.1) 

A(0, x) = eo{x), 5(0, x) = bo{x). 


with initial data / = (eo, 6 o) G (L^(fl))® = 5. Here ^{x) denotes the unit outward 
normal to dfl at a: G F pointing into (,) denotes the scalar product in C^, Utan '■= 
u — {u, v) v, a nd 7 ( 3 :) G C'°°(F) satisfies 7 (a;) > 0 for all a: G F. The solution of the 
problem (1.1) is given by a contraction semigroup (5,5) = V{t)f = t > 0, 

where the generator Gt has domain D{Gb) that is the closure in the graph norm of 
functions u = {v,w) G x satisfying the boundary condition 


Vtan - l{v A Wtan) = 0 On F. 

In an earlier paper [2] we proved that the spectrum of G;, in Re 2 ; < 0 consists of 
isolated eigenvalues with finite m ultip licity. If Gbf = A/ with Re A < 0, the solution 
u{t,x) = V{t)f = e^*f{x) of ( |l.l[ ) has exponentially decreasing global energy. 
Such solutions are called asymptotically disappearing and they are invisible for 
inverse scattering problems. It was proved [5] that if there is at least one eigenvalue 
A of Gb with Re A < 0, then the wave operators W± are not complete, that is 
Ran W- ^ Ran W+. Hence we cannot define the scattering operator S related to 
the Cauchy problem for the Maxwell system and (1.1) by the product W^^W-. 
For the perfect conductor boundary conditions for Maxwell’s equations, the energy 
is conserved in time and the unperturbed and perturbed problems are associated 
to unitary groups. The corresponding scattering operator S{z) : (5^(S^))^ —)■ 
(5^(S^))^ satisfies the identity 


5 -^( 2 ) = S'*(z), zGC 


( 1 . 2 ) 
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Figure 1. Eigenvalues of Gb 


if S{z) is invertible at z. The scattering operator S{z) defined in [S] is such that 
S{z) and S*{z) are analytic in the ’’physical” half plane {z G C : Im^; < 0} and the 
above relation for conservative boundary conditions implies that S{z) is invertible 
for Imz > 0. For dissipative boundary conditions the relation (1.2) in general is 
not true and S{zo) may have a non trivial kernel for some zo,Imzo > 0. Lax and 
Phillips [S] proved that this implies that izo is an eigenvalue of Gb- The analysis 
of the location of the eigenvalues of Gb is important for the location of the points 
where the kernel of S{z) is not trivial. 


The main result of this paper is the following (see Figure 1) 


Theorem 1.1. Assume that for all a: G T, 7 ( 0 ;) ^ 1. Then for every 0 < e 1 
and every N there are constants Cg > 0 and Cn > 0 such that the eigenvalues 
of Gb lie in the region U TZn, where 

A, = {z G C : I Rez| < G,{\lmz\^/^+^ + 1), Rez < 0}, 




N 


C: 1 

Imz < C'Ar( Rez 

+ 1)- 

-N 

5 

B) = 

X{E, B) 7 ^ 0, then 



XE = 

curl B 

on 

D, 

XB = 

—curl if 

on 

D, 

div^ 

= div B = 0 , 

on 

D, 

Etan 

- 7(u A Btan) = 0 

on 

T. 


, Re^ < 0}. 


(1.3) 


This implies that u := {E,B) satisfies 

Au — X^u = 0, on n. 

The eigenvalues of Gb are symmetric with respect to the real axis, so it is sufficient 
to examine the location of the eigenvalues whose imaginary part is nonnegative. 
The mapping z i-A z^ maps the positive quadrant {z G C : Rez > 0,lmz > 0} 
bijectively to the upper half space. Denote by ^/z the inverse map. The part of the 
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Figure 2 . Contours Zi,Z2,Zz,5 = ll2 - e 


spectral domain {A G C : Re A < 0, Im A > 0} is mapped by A = i-^/z to the upper 
half plane {z G C : Imz > 0}. In {z G C : Imz > 0} introduce the sets 

:= {z G C : Rez = 1, /i'^<Imz<l}, 0 << 1, 0 < 5 < 1/2, 

Z 2 := {z G C : Rez = —1, 0 < Imz < 1}, 

Z^ ■= {z G C : I Rez| < 1, Imz=l}. 

Set A = \y/zlh, z G U ^2 U Z 3 . To study the eigenvalues A, |A| > i?o, it 
is sufficient to consider 0 < /i ^ 1. As z runs over the rectangle in Figure 2, 
with 0 < h 1, X sweeps out the large values in the intersection of left and 
upper half planes. The values of z G ^2 near the lower left hand corner, z = — 1, 
of the rectangle go the spectral values near the negative real axis. The spectral 
analysis near these values in Z 2 for dissipative Maxwell’s equations does not have 
clear analogue with the spectral problems for the wave equation with dissipative 
boundary conditions. In fact, for the wave equation if 0 < "/{x) < 1, Va; G T, 
the eigenvalues of the generator of the corresponding semigroup are located in the 

domain (see Section 3, [5] and 0 ). For Maxwell’s equations the eigenvalues of 

Gb lie in the domain A,. U TZn and for 0 < 7 ( 0 ;) < 1 and ^{x) > 1 we have the same 
location (see Appendix for the case K = {x : ||a;|<l}). 


Equation (1.3) implies that on 12 each eigenfunction u = {E,B) of Gb satisfies 

(1.4) 


y/zE = vCurlR, 

1 


y/z B = — vcurl A , 

1 


and therefore (—A — z)E = (—A — z)B = 0. For eigenfunctions {E, B) ^ 0, we 
derive a pseudodifferential system on the boundary involving Etan = E — {E,v)v 
and Enor = {E, v). A semi-classical analysis shows that for z G ZiU Z 3 this system 
implies that for h small enough we have if|r = 0 which yields E = B = 0. By 
scaling one concludes that the eigenvalues A = ^ of Gb lie in the region Aj U Af, 
where 

AI = {z G C : I argz — 7 r| < Tr/d, |z| > i?o > 0, Rez < 0}. 

The strategy for the analysis of the case z G Zi U Z 3 is similar to that ex¬ 
ploited in [5] and [5]. In these papers the semi-classical Dirichlet-to-Neumann 
map JV{z, h) plays a crucial role and the problem is reduced to the proof that some 
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/i—pseudodifferantial operators is elliptic in a suitable class. For the Maxwell sys¬ 
tem the pseudodifferential equation on the boundary is more complicated. Using 
the equation divi? = 0, yields a pseudodifferential system for Etan and Enor- We 
show that if {E,B) 0 is an eigenfunction of Gb, then \\Enor\\H^{T) is bounded 

by Ch\\Etan\\Hl(r)- The term involving E^or then plays the role of a negligible 
perturbation in the pseudodifferentrial system on the boundary and this reduces 
the analysis to one involving only Etan- The system concerning Etan has a diagonal 
leading term and we may apply the same arguments as those of [5] to conclude that 
Etan = 0 and hence Enor = 0. 


The analysis of the case z ^ Z 2 more difficult since the principal symbol g 
of the pseudodifferential system for Etan need not be elliptic at some points (see 
Section 3). Even where g is elliptic, if | Imz| < it is difficult to estimate the 
norm of the difference Oph{g)Oph{g~^) — I- To show that the eigenvalues of Gb 
lying in A4 are in fact confined to the region TZ^ for every N G N, we analyze the 
real part of the following scalar product in T^(r) 

Q{Eo) ■= Re{{N'{z,h) — y/z^)Eo,Eo)L^(r), Eq '■= E\r- 

We follow the approach in based on a Taylor expansion of Q{Eo) at z = — 1 

and the fact that for z = —1 we have Q{Eq) — 0{h^), WN G N. In the Appendix 
we treat the case when K = {x : |x|<l}isa ball and 7 = const. We prove 
that for 7 = 1 the operator G has no eigenvalues in {Rez < 0}, while for every 
7 G M"*" \ {1} we have infinite number of real eigenvalues. 


2. Pseudodifferential equation on the boundary 


Introduce geodesic normal coordinates (yi, y') G on a neighborhood of a point 
a;o G T as follows. For a point cc, y'{x) is the closest point in F and yi = dist (x, F). 
Define u(x) to be the unit normal in the direction of increasing yi to the surface 
yi = constant through x. Thus u(x) is an extension of the unit normal vector to a 
unit vector field. The boundary F is mapped to yi = 0 and 

X = a(yi,y') = /3(y') + yiv{y')- 


We have 


Moreover, 


d _ , n d A dyj d 

dxk ^ dyi ^ dxk dyj ’ 

J—^ 


'^vk{y') 




dyj 

dxk 


{yi,y') 




k = 1 , 2 , 3 . 


j = 1 , 2 , 3 , and 


'^Vk{x)da,J{x) =dy^{f{a{yi,y')). 


/c=l 


Since ||u(x)|| = 1, (y.d^.v) = 0, j = 1,2,3. 

A straight forward computation yields 

u(x) A YCurlM(x) = \hdnUtan+ {{DxiU,v), v)^ 


ihdnUtan + (grad ft,(?x, u) J - ihgo{utan), 

\ / tan 


X G F, 


tan 
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where 


Dxj = -ihd^., j = 1,2,3, grad/,/ = 

9o{'^tan) — {('Il'/an, — 12 , 3 - 

Setting Enor = {E,i'), from one deduces 


V f\B = -^v A vcurli; = -^DyEtan -^ 

^JZ 1 a/z ^/Zl 


hEfiorj tan) 

/ tan 


where Dy = —ihdy and the boundary condition in (1.3) becomes 


( 


Ey ^yZ]Etan 


lEr, 


r') + i/l5o(li'tan) — 0, X S F. (2.1) 

/ tan 


Next 


grad,,/(a;) I tan = ['^-^DyJ{a{yi,y'))'^ 

and for u = {ui ,^ 2 ,^ 3 ) € , 




fc=l, 2,3 


3 3 


o u dy 

jdivu{a{yi,y')) = {Dy^u{a{yuy')),i^{y'))+'^'^ ^Dy^Uk{a{yuy')) 

1 ; 1 • o 

k —1 j —2 


= D 


{^nor (yi, y ~b ^ ^ ^Vj (^^tan lU i ^ , 


yi 


f=2 


where {u{a{yi,y')),v(y')) := Unorivi^y') and Z depends on the second derivatives 
of yj, j = 2.3. Apply the operator Dy-^ — to div i?(a(?/i, y')) = 0 to find 


t /\^vi)^nor{yi,y ) + y~l -Pyj ( (-Pyi ) 






7(y') 


i=2 


7(y')' 


— h{{Dy^ — -^)Etan, -Z^) + h{Etam Zi), 

7 

where 7 (y') := j{l3{y')). 

Taking the trace yi = 0 and applying the boundary condition (2.1 1 , yields 

3 3 


(^yi + T T qI[ ^,)Dy,Ey,or{0,y') 


j,Ai=2 A;=l 

= h 


dyj dyy 

k 

(grad; 


hEjior 


liy'Y 

^^^{0,y'),z) + hQi[Etan[0,y')), 


( 2 . 2 ) 


with 

||Ql(Atan(0,y'))llL2(R2) < C'2||Atan(0,y')llL/i(R2)- 
Here i7^(r), s G M, denotes the semi-classical Sobolev spaces with norm || {hdxYu\\i,- 2 (Y^, 
(hdx) = (1 -I- \\hdx\Yy^^■ In the exposition below we use the spaces (A^(r))^ and 
{Hf^{T)Y of vector-valued functions but we will omit this in the notations writing 
simply L 2 (r) and iF^(r). 


The operator —h'^Ax — z in the coordinates {yi,y') has the form 
V{z, h) = Dl^ E r[y, Dy,) + 91 (y, Dy) + h^q-z 
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with r{y,T]') = {R{y)T]',rj'), qi{y,r]) = {qi{y),r])). Here 




k^l 




is a symmetric (2 x 2) matrix and r{Q,y' ,7]') = rQ{y',ri'), where rQ^y'^rj') is the 
principal symbol of the Laplace-Beltrami operator —h^Ar on T equipped with the 
Riemannian metric induced by the Euclidean one in We have 


(r{z,h)Er,or)iO,y') = {V{z,h)E,v){Q,y') + hQ 2 {E{Q,y')), 


where 


I1Q2(-E(0, j/'))||l2(k 2) < (72 ||E(0, j/')||^i(R2). 

Since V{z,h)E = 0, this lets us replace the terms with all second derivatives of 
Enor in (2.4) by zEnoriQ,y') modulo terms having a factor h and containing first 
order derivatives of Enor- This follows from the for m o f the matrix R{y) given 
above. After a multiplication by the equation (2.2) yields 

(Dy, - -f{y')Vz)EnoriO, y') = hQ 3 (A( 0 , y')), (2.3) 

where Q-i{E{Q,y')) has the same properties as Q2{E{0,y')). 


Let 'ijj{x) G (7“(IR^) be a cut-off function with support in small neighborhood of 
Xq € r. Replace E,B by E^ = Eif;, = Btp. The above analysis works for E-ip 
and B,^ with lower order terms depending on tp. We obtain 

{{D„ - J{x)^/z)E\^^^p{x),v{x)) = hQ 3 ^.^{E\r). 


Taking a partition of unity in a neighborhood of T, yields 

((-Dy - 7(a;)v^)£'|r, v) = hQ4{E\r), ||Q4(L^|r)||L2(r) < (7||E|r||//i(r)- (2.4) 

For z G ZiLI Z 2 U Z 3 let p{x', , z) = \J z — 7’o(a;',^') € C°°{T*T) be the root of 

the equation 

p'^ +ro{x',^')-z = 0 


with lmp{x',^',z) > 0. For large |^'|, 

~ IC'I, 


lmp{x',^',z) ~ l^'l, 


while for bounded |^'|, 

linp{x',^',z) > 

We recall some basic facts about h-pseudodifferential operators that the reader 
can find in [3]. Let X be a (7°° smooth compact manifold without boundary with 
dimension d > 2. Let (x,^) be the coordinates in T*{X) and let a{x,^,h) G 
C°°(T*(X)). Given m G R, I G U, S > 0 and a function c(h) > 0, one denotes by 
Sg‘^{c{h)) the set of symbols so that 

\d^d^^a{x,^,h)\ < C„,^(c(h))-'-^(l“l+l^l)(l + |e|)™-l''l, Va,V/3, ix,0 G T*{X). 

If c{h) = h, we denote S'^’™(c(/i)) simply by Sg^. Symbols restricted to a domain 
where |^| < C will be denoted by a G Sg{c{h)). The h—pseudodifferential operator 
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with symbol a(x, h) acts by 

{Oph{a)f){x) := (27rh)"‘^+^ [ e~^^^~y’^'^/^a{x,^,h)f{y)dyd^. 

Jt*x 

For matrix valued symbols we use the same definition. This means that every 
element of a matrix symbol is in the class Sg'"^{c{h)). 

Now suppose that a{x, h) satisfies the estimates 

\d:aix,^,h)\ < co(h)h-l“l/2, (cr,a e T*(X) (2.5) 

for |a| < d — 1, where co(h) > 0 is a parameter. Then there exists a constant C > 0 
independent of h such that 

\\Oph{a)\\L^(x)^L'^{x) < Cco{h). (2.6) 

For 0 < (5 < 1/2 products of h-pseudodifferential operators are well behaved. If 
a G , b G and s G K, then 

\\Oph{a)Ophib) - Op/i(a5)|l^s(x)^_y=—i-™ 2 +i(x) < ^ 2 . 7 ) 

Let M G be the solution of the Dirichlet problem 

{—h^A — z)u = 0 on n, u = F on F. (2.8) 


Introduce the semi-classical Dirichlet-to-Neumann map 
N{zM : i7/(F) 9 F ^ D^v G 

G. Vodev [9] established for bounded domains iF C d > 2, with (7°° boundary 
the following approximation of the interior Dirichlet-to-Neumann map Afintiz,h) 
related to (2.8), where the equation [—h?A — z)u = 0 is satisfied in K. 

Theorem 2.1 ([5]). For every 0 < e <C 1 there exists 0 < ho{e) <C 1 such that for 
z G Zi^ := {z G Zi, I Imz| > h^~'^} and 0 < h < ho(e) have 

Ch 

\\^^^nt{z,h){F)-Oph{p + hb)F\\Hl^T) < /IT I ll-P^IU^cr), (2.9) 

\J\\mz\ 


where b G Sq ]^(r) does not depend on h and z. Moreover, (2.9) holds for z G Z 2 UZ 3 
with I Imz| replaced by 1. 


With small modifications (2.9) holds for the Dirichlet-to-Neumann map JV{z, h) 
related to ( |2.8| ) (see 0 )- Applying ( |2.9[ ) with A/'( 2 :, h) and F = Eq = E\r, we obtain 

Ch 

< H^olk^y)- 


{Af{z, h)Eo, v) - {Oph{p)Eo, v) 


L2(r) 


Therefore (2.4) yields 


{Oph{p) - 7v^)Fo, v ) - hQi{Eo) 


< 


^/\hr^z\ 
Ch 


L^T) ^llm^l 


\Eo\\L^(r)- 


( 2 . 10 ) 


( 2 . 11 ) 


The commutator [Oph{p),n(x)] is a pseudodifferential operator with symbol in 

^1-5^0.0 

\\[Oph{p),nk{x)]Enor\\Hl(r) ^ C 2 h^~^\\Enor\\Hl(r)’ k = 1 , 2 , 3 , j = 0 ,l. 

h 


The last estimate combined with (|2.11|) implies 

<G3 


{Oph{.p) — l\fz)Enor — hQ^{Ei^) < Gsf' I ll^ollL2(r)- (2.12) 

L^(r) ^, /I Im 7.\ / 


\/| Im; 
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3. Eigenvalues-free regions 

For z G g we have p G 5'°’^ with 0 < (5 = 1/2 —e < 1/2, while for z G Z 2 UZ 3 we 
have p G 5'g’^ (see i)- Since F is connected one has either 7 ( 31 ) > 1 or 0 < 7 (z) < 1. 
We present the analysis in the case where 0 < 7 ( 3 :) < 1, Vx G F. The case 1 < 7 (x) 
is reduced to this case at the end of the section. Clearly, there exists eg > 0 such 
that 

Co < 7 ( 2 ;) < 1 - eo, Vx G F. 

Combing (2.4) and (|2.9[), yields 


ill^o||L2(r) + Ci/i||EollHi(r)j 


\{{Oph{p) - 7 (x)v^)Eo,u(x))||L 2 (r) < C 


\/| Im; 


where for z G .Z 2 U Z 3 we can replace | Imz| by 1. This estimate for Eq and the 
estimate for the commutator [Oph{p),Vk{x)] imply 

I ll^ollL2(r) + (^ 4 / 1 ^ '^ll^ollHi(r)- (3-1) 


\{Oph{p)-l{x)^)Enor\\L^(T)< - i ^ 5| 


\/|Imz| 

Let (x',5') be coordinates on T*(F). Consider the symbol 

c(x',^',z):= p(x',^',z) - 7 (x')v^, x'G F. 

Following the analysis in Section 3, [5], we know that c is elliptic in the case 0 < 
7 {x') < 1 and if z G Zi we have c G | Imz|c“^ G while if z G .^2 U Z 3 

one gets c G S'g’^, c~^ G 5'q’~^. This implies 


\\Oph{c )Ophic)Enor\\Hl(r) < 


C 


Imz 


I ^Ph {,^')^nor\\L'^(V) • 


On the other hand, according to Section 7 in [3], the symbol of the operator 
Oph{c~^)Oph{c) - J is given by 


N 


")ix'+bNix',e) 


3=1 


l“l=i 


:= bN{x',C) + bNix',C), 


where 


|9“,6Ar(x',C')l < 

Taking into account the estimates for c~^ and c, and applying (2.5), and (2.6) yields 
(Oph{c~^)Ophic) - I^E„ 


HUT) 


<C5 


Imz ' 




J = 0,1. 


Repeating the argument in Section 3 in [ 8 ] concerning the case 0 < 7 (x') < 1, for 
z G and 0 < ^ < 1/2, one finds 


\Enor\\HUT) — 


(Oph{c ^)Oph{c) - I^E„ 


]e 

1 ^nor 

+ 

J 

^i(r) 


Oph{c uOph{c)En 


< C6/ll-2^||F;o|U2(r) + C^h^-^^\\Er.,r\\HUT) + || Eg |i (r) ■ 


Clearly, 


HUT) 

(3.2) 


|-^o||H/(r) < W^tanWHUT) + 11 ^"or-11 4/^ (p) ; fc G N 
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with Bk independent of h. Hence we can absorb the terms involving the norms of 
Enor in the right hand side of (3.2 1 choosing h small enough, and we get 

\\Enor\\Hl(T) < ‘^^\\Etan\\Hl{T)- (3-3) 

The analysis of the case z G Z 2 U Z 3 is simpler since in the estimates above we have 
no coefficient | Imz|“^ and we obtain the same result with a factor h on the right 
hand side of (|3.3|). 


With a similar argument it is easy to show that 

\\Enor\\L'^{T) < C \\Etan\\L'^ {T) ■ 


In fact from (2.12) one obtains 
Oph{c~^) ( Oph(p)-7Vz)Enor-hQ4(Eo) 
and 


< 


C. 


L2(r) |Imz| Vy|lmz 


(: 




1-S\ 


(3.4) 


E, 


\\Oph[c ^)QA{Eo)\\L 2 (^y) < 


Co 


Im z| 


iCo||L2(r). 


Combining these estimates with the estimate of ||Op;i(c ^)Oph{c) — I\\L^(r)^L^{r) 
yields (| 0 |). 


Going back to the equation (2.1), we have 


(d^ - -y/^E = (d^ - ^y/^EnorV “ “ l)'/zEnory 

\ ^ / \ / 7 

+ihgoiEtan) + (grad hiEnor)^ a; G T. 


(3.5) 


Notice th at f or the first term on the right hand side of (3.5) we can apply the 
equality (2.4), while for Enor and ^grad/j(C„or)^ we have a control by the 

estimate p.3[). Consequently, setting Eq = E\r, the right hand side of (3.5) is 


bounded by Ch 


1-2(5| 


Co||_f/i(r)- Next 
1 < ^ < 


< —, VxGT. 
eo 


1 - eo 7 (x) 

This corresponds to the case {B) examined in Section 4 of [5]. The approximation 
of the operator N{z, h) given by (2.9) yields the estimate 

\\{Oph{p) -V^)Co||L 2 (r) < c( —=:||Co||L 2 (r) + ^‘^l|Co||/Li(r))■ (3-6) 

For z G U Z 3 the symbol 

(i(x', z) := p(x', i', z) - 

7(x') 


is elliptic (see Section 4, [S]) and d G d ^ G S'^’ ^. Then from (3.6) we esti¬ 
mate ||£'o|l//i(r) and we obtain Eq = 0 for h small enough. This implies E = B = 0. 


Now recall that we have 


Im^z Re^z 

Re A =-,-, Im A = —,—. 

h h 


Suppose that z & Zi. Then 
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ReA| > C(h-^Y-\ |ImA| < Cih-^ < CalReAlT^. 


So if 


ReA| > CallmAI^-^ ReA < -Ca < 0, 


there are no eigenvalues A = of ^b- In the same way we handle the case z € 

and we conclude that if z e U Z 3 for every e > 0 the eigenvalues A = of G{, 
lie in the domain Aj U , where 

A1 = {z G C : I argz — 7 r| < Tr/d, |z| > i?o > 0, Rez < 0}, 

Aj being the domain introduced in Theorem 1.1. Of course, if we consider the 
domain 

^3.5o = {-2 G C : |Rez| < 1, Imz = do > 0}, 
instead of Z 3 , we obtain an eigenvalue-free region with M. replaced by 

Msg = {z G C : I argz — 7r| < arctg do, \z\ > i?o(do) >0, Rez < 0}. 

The investigation of the case z G .^2 is more complicated since the symbol d may 
vanish for Imz = 0 and {xq,Yo) G T*(r) satisfying the equation 


vATAwds) - 


= 0 . 


To cover this case and to prove that the eigenvalues A = with z € Z 2 
are confined in the domain TZn, VA G N, we follow the arguments in and [S]. 
For z G ^2 we introduce an operator T{z,h) that yields a better approximation 
of J\f{z,h). In fact, T{z,h) is defined by the construction of the semi-classical 
parametrix in Section 3, for the problem (2.8) with F = Eq. We refer to [5] for 
the precise definition of T{z,h) and more details. For our exposition we need the 
next proposition. Since (A — z)E = 0, as in (2, we obtain 

Proposition 3.1. For z G Z 2 and every N G N we have the estimate 

\\JViz,h)Eo-Tiz,h)Eo\\Hl^r) < CNh-^°h^WEohnr) (3.7) 

with constants C]\[,so > 0, independent of EQ,h and z, and sq independent of N. 
Proof of Theorem 1.1 in the case z G Z 2 . Consider the system 


Di, — -\fz 1 Etan — 


div hEtan + div h {EnorV^ =0, a: G F, 


^norf ihgoi^E lari'} — 0 ; 

/ tan 


a; G F, 


(3.8) 


where div hF = Yfk=i Ex^Fk- 

Take the scalar product (,)L 2 (r) in T^(r) of the first equation of (3.8) and Etan- 
Applying Green formula, it easy to see that 

Re(grad,7^ian)L^(r) ~ hEtam (3-9) 

tan 

We claim that 

lm{go{E 

tan ), Etan)L^{T) = 0. 

Let Etan = iwi,W 2 ,W 3 ). Then 

3 Qi/ X ^ 1 

{go{Etan),Etan)c3 = ^ ^ ^k^Wj = -{Sw,w)c3, 


k,j=l 


k,j=l 


(3.10) 
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where S := with F(a;) = q{x)iy{x), q{x) > 0 because = 

0. Thus if the boundary is given locally by X 3 = G(xi,X 2 ), we choose V(x) = 
{—dx^G, —dx^G, 1) and it is obvious that S is symmetric. Therefore Im(S'r(;, u>)c 3 = 


0 and this proves the claim. Hence (3.10) implies 


Re[l/l(^Q i^Etan) ^ E tan )L 2 (r)] = 0 . 


(3.11) 


From the -^^(F) scalar product of the second equation in (3.8) with we obtain 

R,e(div, F/yjor)_L^ (r) “1“ ^nor} l^(t) —0- (3.12) 

In fact, 

biv hi^^nor^^ — ^v^nor \hE^Qj.^\Y V 

and Im^div 1= 0 . 

Taking together ( |3.9[ ), ( |3.11[ ) and ( |3.12[ ), we conclude that 


Re 


{{Elf ^ )Eian: Elan) L^(r) “f {EifEnQ^l^-f Enor^) L^(T) 


— HelOiyE, e\ —^G{^^^Etan,Etan)L^(r) — 0. 
\ /L2(r) 7 ^ ' 

Here we have used the fact that 

{DyEianfEjiQxlEjQ^ = Difi^Eiani EjiQxEjQ^ = 0 . 


Applying Proposition 3.1 with i?|r = Eq, yields 
Re(r(z, h)EQ, Eq 


Rs/ Eidu^ Ef^dji 

L2(r) \ 7 / L^{r) 


<GNh-^°h^\\Eo\\LHr)- 

(3.13) 


For z = —1, as in Lemma 3.9 in and Lemma 4.1 in [5], we have 
|Re(T(-l,h)Ao,F;o)L2(r)| < Gwh-«“+'^||F;o||i3(r) = 0. 
Consequently, by using Taylor formula for the real-valued function 


Re 


T(z, h)EQ, Eq 


L2(r) 


Etani Etan , 

7 / L^r) 


we get for every N gN the estimate 
Im 


{E^Zt,h))EQ,EQ 


/ 7i 


i2(r) \2^t 


Etan ; Etan 


L—-sq+AT 


L2(r)J 

2 

L2(r) 


(3.14) 


where = —1 -|- it Imz, 0 < t < 1 , 7 i = 7 

According to Lemma 3.9 in [S], in ( |3.14| we can replace ^{zt, h) by Oph{^^{zt)) 
and this yields an error term bounded by G/i|| . On the other hand, 




Gph{ pf {zt^)Etam Eaor^\ “t”/Op/j( {zt)^ E^qx ^ Etan^\ 

OZ / L^(T) \ OZ / 

|2 


LHT) 


— C'/i||Ao||^2(r) 

since the estimate (3.4) holds for z G Z 2 with factor h and ^{zt) G 
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J := 


Thus the problem is reduced to a lower bound of 

9P, 


Im ( (^Oph ( ^tant ^tan) L^(T)~^{^Ph ( ))-£^nor^i ^nor^) (F) 


> 


lm{(Oph{^{zt)) - Etan, Etan) L'^ir) 


2v^V 


-Ci\\E, 


nor II L 2 (r). 


Since 71 (x) > 1, Vx G T, applying the analysis of Section 4 in [5] for the scalar 
product involving Etan, one deduces 

lui(^(Ophi^izt)) - > 0 - 

By using once more the estimate (3.4), for h small enough we obtain 

J > bl(ll^tan||i 2 (r) + ll-®nor||i 2 (r)^ — Boh\\Etan\\‘L 2 (r) — ^2 |i ^^0 || ^2 (r) i 0 < r ]2 < rji. 

Consequently, (3.14) yields 

F —50+-^ 

(772 — Bih)\\Eo\\j^ 2 (^p'j < Cn I 7 ^ ,,^ll-^olli=’(r) 

and for small h we conclude that for z € Z 2 the eigenvalues A = of Gb lie in 
the region TZn- This completes the analysis of the case 0 < 7 (x) < 1, Vx G T. 


To study the case 7 (x) > 1, Vx G T, we write the boundary condition in (1.1) as 
1 


7(x) 


(u A Elan) (u A (u A Etan')') — / \ ^ Eian') T .Btan — 0. 


7(x) 


Next 


vAE= ^uA vCurlB = -^D^Btan 
yjz 1 Jz 


(grad I 


{Etan) 


Bfar), 


Tgrad hBnor 1 “1“ — 0, X € r (3.15) 

\ / tan 


and one obtains 

(Di. - 7(a^)v^) 

which is the same as (13 with Eban^ Ejiqj- replaced respectively by B^q^ 

and replaced by 7 (x) > 1. We apply the operator Dy^ — 7 -\/z to the equation 
div B = 0 and repeat without any change the above analysis concerning Etam Enor- 
Thus the proof of Theorem 1.1 is complete. □ 

Remark 3.2. The result of Theorem 1.1 holds for obstacles K = where 

Kj, j = 1,..., J are open connected domains with C°° boundary and Ki n Kj = 
0, i ^ j. Let Tj = dKj, j = 1,..., J. In this case we may have 7 (x) < 1 for some 
obstacles Tj and 7 (x) > 1 for other ones. The proof extends with only minor 
modifications. The construction of the semi-classical parametrix in is local and 
for the Dirichlet-to-Neumann map JVj{z, h) related to Tj we get the estimate 


\\Mj{z, h){E) — Oph{p + hb)F\\fji^Y'j) < 


Ch 






The boundary condition in (1.1) is local and we can reduce the analysis to a fixed 
obstacle Kj. If {E, B) 7^ 0 is an eigenfunction of Gb, our argument implies Etan = 0 
for X G Tj if 0 < 7 (x) < 1 on T^- and Btan = 0 for x G Tj in the case 7 (x) > 1 on 
r, . By the boundary condition we get Etan = 0 on T and this yields E = B = 0 
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since the Maxwell system with boundary condition Etan = 0 has no eigenvalues in 
{z S C : Re z < 0}. 


4. Appendix 


In this Appendix, assume that 7 > 0 is constant. Our purpose is to study the 
eigenvalues of Gb in case the obstacle is equal to the ball B 3 = {x £ : \x\ < 1 }. 
Setting A = Im^ > 0, an eigenfunction {E, B) ^ 0 ot Gb satisfies 


cml E =—ifj.B, cuTlB = i^E. (4.1) 

Replacing B hy H =-B yields for {E,H) £ iH^i\x\ < 1))®, 

{ curlA = i^i7, cut\H =—i^E, for x £ B 3 , 

Etan + A Htan) = 0, for X £ S^. 

Expand E{x),H{x) in the spherical functions Y^{uj),n = 0,1, 2,..., \m\ < n, uj £ 
and the modified Hankel functions hn\z) of first kind. An application of Theorem 
2.50 in [3] (in the notation of [3] it is necessary to replace w by ^ S C \ {0}) says 
that the solution of the system (4.2) for |a;| = r = 1 has the form 


Etaniy^') — 


n—1 \m\<n 






1 , 


!•=! |m|<n 


Here U^{uj) = 


=grads 2 F^(a;) and V^{oj) = v A U^(uj) for n S N, —n < 


y'n(n+l) 

m < n form a complete orthonormal basis in 

(S^) = {u £ (L^(S^))^ : (u, u) = 0 on S'^}. 

To find a representation of u A Eltan, observe that u A (u A U^) = —C/™, so 


(uAi?tan)(w) =-7-^ ^ [, 

n—1 \m\<n 

and the boundary condition in (|4.2|) is satisfied if 




ifi) 


= 0, Vn € N, \m\ < n, 


(4.3) 


~ i)r = 0, Vn e N, \m\ < n. (4.4) 

For 7 = 1 , there are no eigenvalues. 

Proposition 4.1. For 7=1 the operator Gb has no eigenvalues in {Rez < 0}. 
Proof. The functions hn\z) have the form (see for example [7]) 

7 «(x) = ,r, 

X ^ m\{2x)'^ in — m)\ x \2x) 

m=0 V / V > 


R«(z) 


z"* (n + m)! 

ml {n — m)\ ™ 

m=n '■ ' m=n 


with 
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Therefore the term in the brackets [...] in (4.3) becomes 

. n , 


m—0 


Setting w = we must study for Re w > 0 the roots of the equation 


1 — 'T 

gn{w) := — — -Rn{w)+ wR'^{w) = 0. 
2w 


(4.5) 


For 7 = 1 one obtains R'^iw) = 0. A result of Macdonald says that the zeros of 
the function hn\z) lie in the half plane Imz < 0 (see Theorem 8.2 in 0 ), hence 
Rniw) 7 ^ 0 for Re w > 0. By the theorem of Gauss-Lucas we deduce that the roots of 
R'niw) = 0 lie in the convex hull of the set of the roots of i?„(w) = 0 , so R'n{w) 7 ^ 0 


for Reu> > 0. Consequently, (4.3) and (4.4) are satisfied only for a™ = = 0 and 

Etan = 0. This implies E = H = 0. □ 


For the case 7 7 ^: 1, there are an infinite number of real eigenvalues. 


Proposition 4.2. Assume that 7 € K“'"\{1} is a constant. Then Gb has an infinite 
number of real eigenvalues. Let 70 = max{ 7 , ^}. Then all real eigenvalues A with 
exception of the eigenvalue 


(7o-i)(i + y^i + ^) 

satisfy the estimate 


max{(7o - 1), V7o - 1} ' 


(4.6) 


(4.7) 


Proof. Assume first that 7 > 1. Then qn{w) = wgn{vu) = 0 has at least one real 
root Wo > 0. Indeed, gn(0) = < 0, qniw) —>■ +00 as w —)■ +00. Choosing 

a^° 7 ^ 0 for an integer toq, |mo| < n and taking all other coefficients a™, fiff equal 
to 0, yields Etan ^ 0 and Gb has an eigenfunction with eigenvalue A = < 0. 

It is not excluded that gn{w) and gm(w) for n 7 ^ m have the same real positive 
root. If we assume that for Rew > 0 the sequence of functions {gn{w)}i^^i has 
only a finite number of real roots wi,..., wjv, Wj € then there exists an infinite 
number of functions gn^ (w) having the same root which implies that we have an 
eigenvalue of Gb with inhnite multiplicity. This is a contradiction, and the number 
of real eigenvalues of Gb is infinite. 

It remains to establish the bound on the real eigenvalues. First, consider the 
case n = 1. Then one obtains the equation 


2w^ 7 ~ 1 

2w + I 2 


which has a positive root wq = 7 ^7 — 1 + 1/(7 — I)^ + 4(7 — I)^. This yields the 
Ai from (14.61) 


Next examine the case n > 2. For a root wq G one has 
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Case 1. Wo > Then the inequality 

Sm=2 ^ + «1^0 


I]m=2 amW™ + OlWo + 1 Z]m=2 + “l^O + 1 


> 1 


is satisfied since 02 = |(n + 2)(n + l)n(n— 1) > 12. Consequently, 2wo < 7 — 1 and 
this implies that the eigenvalue A = — satisfies 


A < - 


7-1 


(4.8) 


Case 2 . 0 < wq < Apply the inequality 


E n m— 1 I 

^^2 'rn-amWQ + ai 


> 


m —1 I 

Z2m=2 0-mWQ + Oi 


E ti m—i I I 1 — v-^n m—i , , i 

^^2 + aiWo + 1 Wo 2^rn=2 O-mW^ + OiWo + 1 


> 2 


that is equivalent to 


(1 - Wo)S'o - OiWo 


“t- di ^ 2 


with So = X)m =2 ' This inequality holds because 

n ^ 

{1 - Wo)amw'^~^ - aiWo > (-a 2 -ai)wo, ai = (n + l)n > 2, 


and, 

-02 — 01 = -(n + 2 )(n+ l)n(n — 1 ) — (n + l)n = n(n+ 1 ) -(n + 2 )(n — 1 ) — 1 


> 0 . 


Therefore, 


This easily yields 


m— 1 

m=l "^Q-mWo 

^■"0 - “'O y^n m , i 

Z^m,= l “m '^0 + ^ 


A < - 


7-1 


, - . (4.9) 

V(7- 1) 

In the case 0 < 7 < 1 one has I /7 > 1 and we apply the above analysis to the 


equation (4.4). Setting 70 = max{ 7 , 4} and taking into account (4.8) and (4.9), we 


obtain the result. This completes the proof. 


□ 


Remark 4.3. Proposition 4.2 yields a more precise result than that in [I] since 
we prove the existence of an infinite number of real eigenvalues Gj, for every 7 S 
K+ \ { 1 }. In the case 7 = e > 0 the eigenvalue Ai has the form 


Ai = 




and this result for small e > 0 has been obtained in [T]. Clearly, as 7 —>■ 1 the real 
eigenvalues of G;, go to — 00 . . 

It is easy to see that for 7 > 1 the equation 5 „(w) = 0 has no complex roots. 
Denote by 

Zj, Kezj <0, j = I, ...,n, n > I 
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the roots of Rn{w) 
Then 


and 


0. Suppose that gniwo) = 0,n > 1 with Remo > 0, Immo ^ 0. 


Im 


' 1-7 
- 2mo 


n 


i=i 



= 0 


(1 - 7 ) Immo 

2KP 


Re wq 


-E 


Im Wq 


^ |mo - 


Im Zj 
^ |mo -Zj p 


E 


T ^0 — r^ / ^ 1 /A\ 

+ ImmoV I .2 = 0- 4.10 

On the other hand, if Zj with Imzj ^ 0 is a root of i?„(m) = 0, then zj is also a 


root and 


Im Zj 


Im; 


Im Zj 


|mo - . 


\wo - ■ 


Equation ( 4.I0[ ) becomes 

^ 7-1 


Im Wq 


-E 


Re^j 


2kop 


|mo - ^iP|mo - 

4Im wo{lmZj)‘^ 

ko - ^iPko - %P 

E 


(ko - Zjl"^ - |mo - Zjf) 


4Remo(Im2;jk 


|mo - ^iPko - Zj 


= 0 . 


(4.11) 


The term in the brackets [...] is positive, and one concludes that Immo = 0. 

Repeating the argument of the Appendix in [ 8 ] , one can show that for 0 < 7 < 1 
the complex eigenvalues of lie in the region 


|z G C : |arg z — 7 r| > 7 r/ 4 , Rez < o|. 


Remark 4.4. We do not know if there exist non real eigenvalues for B 3 . 

Acknowledgment. We thank Georgi Vodev for many useful discussions and 
remarks concerning an earlier version of the paper. 
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